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H ' ABSTRACT 



We have measured the equivalent width of the Ha emission hne for 11006 galax- 
ies brighter than Mb = -19 (17a = 0.7, 17„ = 0.3, Hq = TOkms-^Mpc"^) at 
0.05 < z < 0.1 in the 2dF Galaxy Redshift Survey (2dFGRS), in the fields of seventeen 
known galaxy clusters. The limited redshift range ensures that our results are insen- 
sitive to aperture bias, and to residuals from night sky emission lines. We use these 
measurements to trace /i*, the star formation rate normalized to L* , as a function of 
distance from the cluster centre, and local projected galaxy density. We find that the 
distribution of /i* steadily skews toward larger values with increasing distance from 
the cluster centre, converging to the field distribution at distances greater than ^ 3 
times the virial radius. A correlation between star formation rate and local projected 
density is also found, which is independent of cluster velocity dispersion and disap- 
pears at projected densities below ~ 1 galaxy (brighter than = —19) per Mpc-^. 
This characteristic scale corresponds approximately to the mean density at the cluster 
virial radius. The same correlation holds for galaxies more than two virial radii from 
the cluster centre. We conclude that environmental influences on galaxy properties are 
not restricted to cluster cores, but are effective in all groups where the density exceeds 
this critical value. The present day abundance of such systems, and the strong evolu- 
tion of this abundance, makes it likely that hierarchical growth of structure plays a 
significant role in decreasing the global average star formation rate. Finally, the low 
star formation rates well beyond the virialised cluster rule out severe physical pro- 
cesses, such as ram pressure stripping of disk gas, as being completely responsible for 
the variations in galaxy properties with environment. 
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1 INTRODUCTION 

The effect of local environment on galaxy evolution in gen- 
eral is not well understood. Studies of environmental effects 
in the past have been largely devoted to the study of galaxies 
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more common galaxies (e.g. 


Dressler 1980; Dressier et al. 


1985; C 


ouch & Sharpies 1987; 


Balogh et al. 1997, 1999; Pog- 


giant i c 


t al. 19991; 


Moss & Whittle 200C; Couch et al. 200ll 


Solanes et al. 2001 


). However, galaxies in cluster cores com- 



prise only a small fraction of the stellar content within the 
universe, and thus it is not obvious that the processes which 
effect these galaxies are important for galaxy evolution in 
general. 

More recently, however, work has begun to show that 
star formation is suppressed in cluster galaxies far from the 
core. From the CNOCl cluster sample, Balogh et al. (1997 



199S) found that the mean cluster galaxy star formation 
rate may be suppressed as far as twice the virial radius (Rv) 
from the cluster centre, relative to a field sample selected in 
the foreground and background of the clusters. However, the 
data at large radii were sparse, and being derived from only 
a few clusters were sensitive to the effects of substructure 
and non-sphericity. Thus it is not possible to draw strong 
conclusions about the relative cluster galaxy star formation 
rate beyond the i?„ from these data. Wide field photometric 
analysis of clusters using Subaru has recently suggested that 
the tight red sequence of early type galaxies first presents 
itself in small groups of galaxies within the infall region of 
the mas sive cluster C10939-)-47 a,t z = 0.39 (Kodama et al 



200l|)7]rhis is the first work to suggest that a "critical" en- 
vironment for galaxy evolution exists. A larger survey, de- 
signed specifically to study the outer regions of clusters is the 
Las Campanas/ Anglo Australian Observatory Rich Cluster 
Survey (LARGS), a sample of 17 rich. X-ray bright clusters, 
with photometry and spectroscopy extending out to very 
large radii (~ 6 Mpc). Early results confirm the radial gra- 
dient in photometric and spectrosco pic properties out to the 
virial radius and, perhaps, beyond ( p'Hely 200C ; Pimbblet 
et al. 2l|)01a|,[b|). 

It therefore seems likely that galaxy star formation rates 
are reduced before they are accreted by a cluster, for exam- 
ple in smaller groups. If this is the case, the imphcations 
could be p rofound, as most ga l axies at the present day are 
in g roups ([Turner fc Gott 1976| ; [Oeller fc Huchra 1983| ; [Tully 
1987; Qarlberg et al. 2001 ); if environmental processes are 
important in these regions, they will clearly be reflected in 
the evolution of the universe as a whole. As structure builds 
up in the universe, more and more galaxies can be found 
in groups and, if these environments serve to terminate star 
formation, the mean star formation rate of the universe will 
decline. This might explain at least part of the o bserved de- 



clin e in global star form a tion with cosmic t ime (Lilly et al 
199e |; [K^adau et al. 199^ ; [Cowie et al. 1999| ). 



The 2dF galaxy redshift survey (2dFGRS) allows the 
unprecedented opportunity to study the spectroscopic prop- 
erties of galaxies at an arbitrarily large distance from any 
given clust er. The details of t he survey strategy are given 
elsewhe re (Colless et al. 2001 ), but summarized briefly in 
Section 2.1, Analysis of the whole sample will allow a defini- 
tive study of any correlation between spectral properties (i.e. 
emission line strength) as a function of a continuous variable 



like local density. For this preliminary study, we are specif- 
ically interested in establishing precisely where galaxies in 
the vicinity of known clusters begin to exhibit properties 
which differ from those of the average galaxy. We base this 
on a sample of 17 kn own rich clusters within the 2dFGRS, 
from the catalogue of De Propris et al. (2002). 



Our cluster selection, galaxy sample, and star formation 
rate measurements are described in Section ^. In Section ^ 
we show the trend of increasing star formation activity with 
both increasing cluster-centric distance, and decreasing local 
projected density. This is compared with numerical mod- 
els in Section ^ We summarize our findings in Section ^ 
Throughout this paper, we use a cosmology with SIa = 0.7, 
fim = 0.3, Ho = 70kms"^Mpc"\ We use the symbol Mt to 
denote absolute magnitudes measured in the 2dFGRS pho- 
tographic blue system. 



2 DATA ANALYSIS 
2.1 Spectroscopic Data 

The 2dF Galaxy Redshift Survey has obtained over 220 000 
spectra of galaxies located in two contiguous declination 
strips, plus 99 randomly located fields. One strip is in 
the southern Galactic hemisphere and covers approximately 
80° X 15° centred close to the SGP. The other strip is in 
the northern Galactic hemisphere and covers 75° x 10° . The 
99 random fields are located over the entire region of the 
APM galaxy catalogue in the southern Galactic hemisphere 
outside of the main survey strip. Full de tails of the survey 
strategy are given in Colless et al. (2001). 

The survey spectra cover the wavelength range 3600- 
8OOOA at 9A resolution. Only the wavelength range of 3600- 
77OOA is used during the line fitting procedure due to poor 
signal to noise and strong sky emission in the red part of the 
spectrum. The wide wavelength range is made possible by 
the use of an atmospheric dispersion compe nsator (ADG) 
within the 2dF instrument (Lewis et al. 2002). 



2.2 Cluster Selection 



We se lect 17 clusters from the catalogue of De Propris et al 
2002|), in which cluste rs from the Abell catalogues 



bell 



195^ ; lAbeU et al. 19891), the AP M ( palton et al. 1997| ) and 
the EDGG QLumsden et al. 199^ ) were cross-referenced with 
the 2dFGRS. This catalogue is still partially incomplete, but 
the completeness is generally better than 75% within ~ 5 
Mpc of the cluster centres. The mean redshift and velocity 
dispersions of the clusters in this catalogue have been recom- 
puted from the 2dFGRS spectra, and the cluster centroid is 
taken to be the brightest cluster galaxy with early-type mor- 
phology, identified from POSS plates. 

For this analysis, we extract from the 2dFGRS all galax- 
ies within ~ 20 Mpc of the centre of 17 clusters, selected to 
lie at 18000 kms"^ < cz < 29 000 kms"^ The lower veloc- 
ity bound is chosen to limit the angular size to a reasonably 
small, manageable value; the upper limit is defined as the 
velocity at which Hq is redshifted into the first set of strong 
night-sky OH emission lines. Ten clusters were selected to 
have velocity dispersions a > 800 km/s, while the remaining 
seven are systems with 400 km s~^< a < 800 km s~^. The 
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redshift histograms for the 17 clusters, including all galaxies 
brighter than Mb = — 19 within 5 Mpc (projected) of the 
centre, are shown in Fig. 0. Details of the clusters, includ- 
ing their redshifts (cz), velocity dispersions (a), number of 
cluster members brighter than Mb = —19, and complete- 
ness (within 5 Mpc) , are summarized in Table |l| [Pe Propris 
et al. ( ?002 ) resolved Abell 1238 into two clusters aligned 



along the line of sight; we here consider the lower redshift 
cluster, designated Abell 1238L. The cluster centres and ve- 
locity dispersions are generally better determined than they 
appear in Fig. |l|, as they are computed including fainter 
galaxies over a smaller projected area (where the contrast 
with the field is greater). 

2.3 Ha measurements 

All of the measurements of equivalent width have been per- 
formed using a completely automatic procedure. For each 
spectrum we remove the continuum by subtracting the me- 
dian over a ISSA (31 pixel) wide window after first excluding 
known absorption and emission line regions by making use 
of the known galaxy redshift. Bad pixels and sky line residu- 
als and the atmospheric and fibre absorption bands are also 
excluded from the continuum fitting. 

Both emission and absorption lines are fitted with Gaus- 
sian profiles which are adequate for most of the emission 
lines and cores of the absorption lines. Up to 20 individual 
absorption and emission lines are fitted simultaneously using 
a modified Levenberg-Marquardt algorithm. The width and 
height of each line are fitted together with a small perturba- 
tion of the observed redshift. Some lines were constrained to 
be emission or absorption. Others were allowed to be either. 
Note that the relative wavelength spacing of all lines remains 
fixed, but the fitted redshift is allowed to vary slightly (typ- 
ically Az ~ 0.00025, and always Az < 0.005). By fitting 
many lines simultaneously we avoid individual line fits shift- 
ing to the nearest available peak or dip in the spectrum. By 
fitting both absorption and emission lines we ensure that 
the method is robust to the redshift solution whatever type 
of spectrum is being fitted. 

With this technique of simultaneous line fitting it is 
possible to allow for line blends by simply requesting two or 
more lines to be fitted to the blend. For example H/3 is best 
fitted by a combination of a narrow emission and a broad 
absorption line, and the Ha emission line can be accurately 
deblended from the adjacent [Nii]A6548A and [Nii]A6583A 
lines, despite the 9A resolution of the spectra. The [Nil] lines 
are constrained to be in emission while the Ha line may be 
either emission or absorption. To fit the Gaussian profile 
to the data points a consideration has to be made for the 
effect of the undersampling of the data. The solution is to 
model a Gaussian profile which, when undersampled, fits the 
observed data closely. Fig. |2| shows the resulting fit for four 
spectra with varying [Nil] /Ha ratios, and demonstrates the 
effect of the undersampling. 

After line fitting, the parameters of the fit (amplitude, 
sigma and area) and the rms residuals are used to classify 
the quality of the line fit. Usual reasons for rejecting a fit 
are if the line is too narrow (e.g. a noise spike or residual 
cosmic ray hit), or too broad (for a forbidden line). Some 
combinations of lines are also rejected, for example Ha ab- 
sorption combined with [Nil] emission. Lines which are too 



weak for a good fit are also fiagged; however they are not 
rejected from the analysis, so that we retain a dispersion in 
the flux which reflects the measurement uncertainties. Par- 
tial failures are fiagged often due to large rms residuals to 
the fit (broad non-thermal emission lines are poorly fitted 
by Gaussian profiles) or a poor wavelength calibration at 
the blue end of the spectrum, which can lead to a poor line 
profile for [On]. The latter is usually the case when the ob- 
served spectrum was close to the edge of the CCD. Care is 
also taken when a bad pixel has been masked out from the 
spectrum within 2a of the line centroid. 

Equivalent widths are then simply calculated using the 
continuum fit and the measured line flux. A small number 
of spectra are degraded by poor sky subtraction at the data 
reduction stage, which can result in a negative continuum, 
making the EW meaningless. These cases can be easily re- 
moved from subsequent analysis. 

2.4 Sample Selection and Star formation Rates 

2dFGRS spectra within 20 Mpc of each cluster centre were 
extracted from the database. The extreme ends of spectra 
taken earlier than A ugust 1999 are sev erely affected by prob- 
lems with the ADC ( Lewis et al. 200^ ; hence we restrict our 
analysis to data taken after this date. This leaves us with 
53018 galaxies, but to limit the effect of aperture bias and 
sky-subtraction residuals, we restrict the sample to those 
galaxies which lie within 0.05 < z < 0.1. Within this red- 
shift range, the galaxy sample is complete to Mb — —19, and 
we adopt this as our luminosity limit. This leaves us with 
12020 galaxies. For computations of star formation rates, 
we exclude galaxies in which the continuum was neg a,tive , 
or a Gaussian was a poor fit to the line (see Section 2.3). 
This removes an additional 734 galaxies from the sample 
(~ 6%). Finally, for galaxies with a significant Ha equiva- 
lent width (W^Hci > 10 A) we exclude galaxies in which the 
equivalent width of the adjacent [Nll]A6583 line is greater 
than 0.55WHa. These 280 galaxies (2.3% of the sa mple) are 
likely to have a sign ificant non-thermal component (Veilleux 
& Osterbrock 1987). This leaves us with a final sample of 
11006 galaxies. We take cluster members to be those within 
3(T, where a is the cluster velocity dispersion determined 
by [De Propris et al. (200^ , shown in Table |l|. The num- 



ber of such members within the virial radius (see below) is 
denoted A^mom in the table. Note that for the three highest 
redshift clusters A0933, A0954 and A1189, the highest ve- 
locity members are not included due to our overall redshift 
cut (0.05 < z < 0.1); however, all galaxies within 2a are still 
available. 5829 galaxies in our final sample are thus deflned 
as cluster members. 



Kennicutt (1983, 1992) derived a conversion from Ha 
luminosity to star formation rate, under the assumptions of 
Case B recombination, no escape of Lyman-a photons, and a 
Salpeter-like initial mass function. This may underestimate 
the current star formation rate by a s mall factor, due to ex- 
tinct ion in the line-emitting regions (Chariot & Longhetti 
2001). Also, if the nature of star formation is burst-like, the 



instantaneous star formation rate may not be represent ative 
of the average o ver even short (~ 100 Myr) timescales ( Suf 



livan et al. 2001). However, neither of these effects are likely 
to affect a comparison of galaxy populations with similar 
luminosity functions, as is the case in the present work. 
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Figure 1. Redshift histograms for the seventeen clusters used in this study. Only galaxies brighter than Mj, = — 19 and within 5 Mpc 
of the cluster centre are included. The dotted line overlayed on each histogram represents a Gaussian with a central redshift (cz) and 
velocity dispersion as tabulated in Table hi Left: The ten clusters with velocity dispersion a > 800 km s~^. Right: The remaining seven 
clusters with a < 800 km s~^. 



Table 1. The seventeen clusters used in this study 



Name 


R.A. 


Dec. 


cz 


^ mem 


cr 


Completeness 


Rv 


Rv (alt.) 




(B1950) 


(km s-i) 




(km s-i) 


(within 5 Mpc) 


(Mpc) 


(Mpc) 


S0258 


02:23:33.21 


-29:50:26.9 


18060 


31 


583 


0.72 


1.6 


1.9 


ED652 


02:25:11.88 


-29:51:00.7 


18001 


21 


564 


0.75 


1.4 


1.8 


A3094 


03:09:16.42 


-27:07:08.4 


20475 


63 


774 


0.84 


2.0 


2.4 


S0333 


03:13:04.34 


-29:25:41.3 


20042 


40 


998 


0.90 


1.6 


3.2 


S0340 


03:17:55.68 


-27:11:45.6 


20281 


18 


939 


0.87 


1.2 


3.0 


A0933 


10:05:14.50 


+00:45:25.7 


29180 


72 


420 


0.54 


2.4 


1.3 


A0954 


10:11:11.10 


+00:07:40.2 


28622 


74 


832 


0.77 


2.2 


2.5 


A1189 


11:08:30.14 


+01:21:42.6 


28824 


51 


814 


0.77 


1.9 


2.5 


A1200 


11:10:03.25 


-02:56:27.6 


24970 


38 


825 


0.83 


1.7 


2.6 


A1238L 


11:20:20.36 


+01:23:19.4 


22160 


53 


586 


0.82 


1.9 


1.8 


A1620 


12:47:29.78 


-01:16:07.1 


25513 


51 


1095 


0.89 


1.8 


3.4 


A1651 


12:56:47.48 


-03:55:36.9 


25152 


46 


817 


0.47 


2.1 


2.5 


A1663 


13:00:18.05 


-02:14:57.7 


24827 


75 


884 


0.80 


2.1 


2.7 


A1692 


13:09:41.25 


-00:39:59.7 


25235 


49 


686 


0.80 


1.8 


2.1 


A 1750 


13:28:36.52 


-01:28:15.9 


25647 


83 


981 


0.62 


2.4 


3.0 


ED119 


22:13:32.57 


-25:55:10.7 


25546 


38 


1112 


0.84 


1.7 


3.4 


S1086 


23:02:06.51 


-32:49:14.8 


25605 


74 


502 


0.53 


2.4 


1.5 



Since the 2dFGRS spectra are not flux calibrated, we 
cannot derive Hq luminosities, or star formation rates. How- 
ever, after making a small (2A) correction for the underlying 
stellar absorption, we can use the equivalent widths to cal- 
culate the star formation rate normalised to a fiducial lumi- 
nosity (essentially a star formation rate per unit normalized 
luminosity). If p is the star formation rate in units of Mq 
yr ~^ and Z/Hq is the total luminosity of the Ha emission 
fine in ergs s^^, we can define 



ri = fi/Luo 



(1) 



We will use the "average " conversion fac tor of r; = 7.9 x 
10"*^ Mq s yr"^ ergs"^ ([Kennicutt 199^). The equivalent 



width of Hof, corrected for stellar absorption, is given by 

Whc ~ iHc/Lc, (2) 

where Lc is the continuum luminosity in units of 
ergs s~^A~^. We can then calculate fi* as 



LJL* 



T]WHaL*, 



(3) 



where L* is a characteristic luminosity, for normalisation, 
in units of ergs s^^A"'^. We take L* to correspond to the 
knee in the luminosity functi on in the r' band (n ear rest- 
frame Ha), as determined by Blanton ct al. (2001), Mr = 
-21.8 ( Ha = 0.7, f7„ = 0.3, h = 0.7), or L* = 1.1 x 10^" 
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Figure 2. Four examples of the Gaussian line-fits to the spectra, for varying strengths of [Nil] and signal-to-noise ratio. Each plot shows 
the observed low resolution data as a series of filled circles with rms error bars. The modelled Gaussian fit is shown as a smooth curve 
and the modelled Gaussian sampled at the same resolution as the observed data is shown as a histogram. 



ergs s ^ . Therefore, we have 



Integr ating the best fit Schechter function from Cross et al 



^ 

which gives the star formation rate, in units of MqYT^^ , 
normalized to L*. 

We measure the projected distance of each galaxy from 
the cluster, as defined by the brightest central galaxy. In 
some cases the cluster membership of a galaxy is ambiguous, 
because it lies within 20 Mpc and the 3cr redshift limits of 
more than one cluster (e.g. clusters S0258 and ED652). In 
this case, the galaxy is assumed to belong to the cluster 
which is nearest in projected distance. 

In order to put all the clusters (which span more than a 
factor of two in velocity dispersion) on a common scale, and 
to facilitate comparison with theory, we need to relate pro- 
jected distances to the virial radius, _R„, of the cluster. We 
show the spatial distribution of the cluster members within 
1 degree of the centre for each cluster in Fig. |^. From this fig- 
ure it is evident that many of our clusters are not spherically 
symmetric. Thus we must be cautious in our interpretation 
of Rv as a physically meaningful scale, particularly when 
considering individual clusters. 

The definition of i?„ is 



p(< Rv) = Ac(z)pc(z) = Ac(z)pt{z)/il„^{z), 



(5) 



where p(< Rv) is the mean cluster mass density within 
Rtj, Pc and pb are the critical density and mean back- 
ground mass density, respectively, and Ac is the redshift- 
dependent contrast parameter, determined from spherical 
collapse theory. For a flat Sim ~ 1 universe, Ac = 178; 
for our ado pted cosmology at z = 0.07, Ac ~ 107 (Eke 
et al. Ip96|), and Qm{z) = 0.343, so Ac{z)/D,m{z) = 312. 



We will assume that the number density of galaxies is di- 
rectly proportional to the dark matter density, independent 
of scale or galaxy luminosity. In this case, we can take the 
mean background density pt from the luminosity function. 



(2001), we find that the number density of galaxies brighter 
than Mb = -19 is pb = 0.0076 Mpc"^ {h = 0.7). We de- 
termine p(< Rv) by counting the number of cluster mem- 
bers A'' within (weighting by the completeness given 
in Table ^ and assuming a spherical cluster geometry, so 
p(< Rv) = 3A^/(47r_R^). Substituting this into Equation ||, 
we need to solve Rv — 0.465A^'^^. This is done iteratively, 
by first estimating Rv, counting the number of members N 
within Rv , and then recomputing Rv . This is repeated until 
the solution converges, usually within ~ 3 iterations. These 
measurements of Rv are given in column 7 of Table ^. 

Alternatively, the virial radius can be determined di- 
rectly from the velo city dispersion, und er various assump- 
tions, as outlined in Girardi et al. (1998). If M„ is the viri- 



alised mass, and Rv is the cluster virial radius, we have 

3Mv 



'i-KpcKi ' 



(6) 



The virial mass can be related to the velocity dispersion 
a and Rv under the assumption of spherical symmetry, 
through 



so we have 



Rv 



Mv = 3G 



3a 
AnGpcAc 



-1 2 D 

cr Rv, 



-rr/Ho. 



(7) 



(8) 



For either a flat, flm = 1 cosmology (with h — 0.5) or the 
f^A-dominated cosmology we have adopted {h = 0.7), the 
virial radius in Mpc is Rv ~ 3.5cr(l -I- z)~^'^, for a in units 
of 1000 kms-\ 

For 9 of the 17 clusters, this calculation (listed as Rv 
(alt.) in Table ^) agrees with the previous one to within 
~ 20%. For most of the remaining cases, where there is 
a large discrepancy between the two measurements of Rv, 
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Figure 3. Spatial distributions for cluster members within 1 degree of the cluster centre. The filled circles are galaxies with Wna > 
20A. The large circle in each panel traces the estimated virial radius for the cluster. Left: The ten clusters with velocity dispersion 
cr > 800 km s~^. Right: The remaining seven clusters with cr < 800 km s~^. 



the velocity histograms are significantly non-Gaussian, and 
thus the velocity dispersion is likely to be a poor tracer of 
the mass. For this reason, we will always adopt the first cal- 
culation of Rv as the most likely to be correct. Moreover, 
this occaisional discrepancy, and the non-Gaussianity of the 
corresponding velocity histograms, likely implies that the 
computed velocity dispersions are not always simply related 
to the virialised mass. For example, some clusters (ED119, 
S0333, S0340) may have velocity dispersions which are artifi- 
cially inflated by the presence of foreground and background 
structures. Thus, our division of the sample into two based 
on velocity dispersion may not reflect a perfect division into 
low- and high-mass clusters. 

We will draw the reference field population from the 
2400 galaxies more than 6cr from the cluster redshift; i.e. in 
the foreground and background of the clusters. Due to the 
small redshift range considered, 0.05 < z < 0.1, and the 
use of an absolute luminosity limit, the field sample is also 
volume limited. The luminosity function of the field sample 
is comparable to that of the cluster sample, as shown in 
Fig. I (see also |De Propris et al. 2002| ). 



3 RESULTS 

3.1 General Cluster Properties 

In Fig. Is] we show the distribution of normalized star for- 
mation rate, /i*, in the cluster and field samples, excluding 
gala xies with relatively strong [Nn]A6583 emission (see Sec- 
tion 2.4). The cluster sample is limited to the 440 members 



within Rv, while the field sample is drawn from the 2400 
galaxies beyond 6a in velocity. The difference between the 
distributions is highly significan^, with the field galaxy pop- 
ulation weighted toward galaxies with stronger star forma- 
tion. 



3.2 Radial Dependences 

It is well known that star formation ac tivity in clusters in- 
creas e s wit h distance from the centre (Balogh et al. 1997, 
1999). In Fig. |^ we show how the mean and median 



199S 



value of /i* depend on radius in our cluster sample, and com- 
pare that with the field value. We also show, in the third 
panel, the fraction of galaxies with fi* > 1, which represent 
the tail of the distribution, comprised of galaxies that are 
currently forming stars at a high rate relative to their lumi- 
nosity. The sample is also broken up into clusters with high 
(cr > 800 km s~^, triangles) and low velocity dispersions 
(cr < 800 km s~^, crosses). The properties of the field sam- 
ple are shown as the horizontal, solid line. The dashed lines 
bracketing the field line represent the 1-cr standard deviation 
from field to field, computed by ordering the field galaxies 
in right ascension and treating every 200 galaxies as an in- 
dividual sample. This gives some estimate of the expected 
cosmological variance in the field value. 



^ The probability that the two distributions are not drawn 
from the same population is > 99.999% as determined by a 
Kolmogorov-Smirnov test. 
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Figure 4. Luminosity functions of the cluster sample (solid his- 
togram) and the field sample (dashed histogram). 



All three statistics demonstrate that the cluster distri- 
bution of /i* becomes equivalent to the field value only well 
outside the virial radius, at i? <; 3/?^, in excellent agree- 
ment wi th preliminary resu lts from the Sloan Digital Sky 
Survey (Glomez et al. 2002). The implications of this are 
that a representative sample of field galaxies cannot be ob- 
tained within <^ 6 Mpc of the cluster core. Thus, photometric 
studies of clusters which attempt a statistical background 
subt raction by taking the field from the cluster o utskirts 



(e.g. Kodama & Bower 2001; Pimbblet et al. 20011:) are not 



subtracting enough star-forming galaxies, and artificially in- 
flating the number of blue galaxies within the cluster. 

Since many of the clusters are not spherically symmet- 
ric, the interpretation of radial gradients, and the physical 
meaning of Rv, is not straightforward. From Fig. ^ it is clear 
that there is often considerable structure, both within and 
without the virial radius. Furthermore, the galaxies with 
strongest Ha emission (solid points in Fig. |3|) appear to 
be spread evenly throughout the field, avoiding the densest 
regions, regardless of clustercentric distance. Thus, in the 
following section we consider the correlation between star 
formation rate and local density. 

3.3 Density Dependences 

There has been controversy over whether or not galaxy pop- 
ulations correlate mo st closely with c luster-centric radius 
( Whitmore et al. 1993| ) or local density ( Dressier 1980| ; 



Post- 



man fc Geller 198^ If radius is the primary determinant 
anywhere, it is most likely only within the very central re- 
gions of the cluster (Dominguez et al. 2001). Studies which 



stack many clusters to approximate a spherically symmetric 
supercluster circumvent this difficulty, since average den- 
sity becomes a moii o tonic function of radius within Ru (e.g. 
Balogh et al. 1997, 1998). In our case, the outer regions 




of the clusters often contain several large groups or other 
clusters of galaxies (see Fig. ra). Thus, it is probably more 



/X* (Mgyr-i) 

Figure 5. The distribution of star formation rate per unit lu- 
minosity, in the cluster and field samples. The cluster sample is 
limited to galaxies within the virial radius. 



appropriate to consider the local density of the galaxies as 
the most physically interesting variable. 

To compute the local density of cluster members, we 
consider all galaxies in the spectroscopic catalogue (includ- 
ing those with bad ADC or Ha measurements) brighter than 
Mi, = —19, and within Scr of the cluster redshift. We then 
take the distance to the tenth nearest galaxy, in projected ra- 
dius, as rio; the local projected density is then E = lO/vrriQ. 
For galaxies near the boundary of a cluster catalogue, this 
will underestimate the true density. To partially account for 
this, we only consider galaxies within 18 Mpc of the cluster 
centre, so they are at least 2 Mpc from the edge of the cata- 
logue. In some cases, however, the current 2dFGRS database 
is incomplete within the 20 Mpc extracted area, and the den- 
sities of galaxies near these incomplete regions will still be 
underestimated. 

In Fig. ^ we show the distribution of density, for galaxies 
in three radial bins. In the cluster centre, almost all galaxies 
are in regions of very high local density. However, at large 
radii galaxies can be found in a wide range of environments; 
in particular it is not uncommon to find galaxies at _R > 3-Ri, 
with local densities as large as those within the virialized re- 
gion. Within the virial radius, the distribution of E is sim- 
ilar for both high and low velocity dispersion clusters; the 
means are the same within ~ 5%, and the probability that 
both distributions are drawn from the same population is 
0.12 as determined by a Kolmogorov-Smirnov test. Between 
1 < iii, < 5, however, there is a significant (> 99.999%) 
difference, and the mean local density of the clusters with 
(T > 800 km s~^ is more than twice as large as that of the 
lower velocity dispersion clusters. 

In Fig. H, we show the properties of the cluster ^* dis- 
tribution, as in Fig. ^, but plotted against E. The vertical 
line shows the mean projected density of galaxies within the 
virial radius, A'^(< R„)/-kR^ (note that this is not the same 
as the average of the E values calculated for each galaxy) . As 
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Figure 6. The mean {bottom panel) and median (middle panel) value of /i* in the cluster sample, as a function of radius. In the top 
panel we show the fraction of galaxies with fi* > IMq yr~^. Error bars are a jackknife resampling estimate. Solid points represent the 
full galaxy sample, while the triangles and crosses represent only the clusters with a greater than or less than 800 km s~^, respectively 
(offset for clarity). Only points in which the radial bin contains at least three galaxies are shown. The horizontal, solid line represents 
the value of each statistic in the field sample. The dotted lines which bracket the line are an estimate of the 1—cr field to field standard 
deviation, for independent samples of 200 galaxies. 



in Fig. ^ the horizontal lines show the values of each statis- 
tic in the field. The field actually spans a range of densities, 
likely similar to that seen far (> from the cluster cen- 

tre (see Fig. however, our density estimate in clusters is a 
measure of the galaxy density projected along a line-of-sight 
column of unknown length, and thus cannot be directly ap- 
plied to field galaxies to obtain a comparable measurement 
of local density. Thus, it is evident that star formation is 
suppressed at densities of S ~ 1.5 galaxies Mpc~^, approx- 
imately 2.5 times lower than the mean projected density of 
the cluster virialized region. 

As in Fig. ^ the star formation rate distribution at a 
given density is similar in both high- and low- velocity disper- 
sion clusters. This suggests that galaxy star formation rates 
depend only on the local density, regardless of the larger- 
scale structure in which they are embedded, although we 
repeat our caution that the velocity dispersions may not be 
directly related to the cluster mass in all cases. Furthermore, 
as we show in Fig. the correlation of star formation rate 
with density holds at r > 2_R„, well outside the virialised 
cluster region. This demonstrates that star formation is low 
relative to the global average in any region exceeding the 
critical density of 1 galaxy (brighter than Mb = —19) per 
Mpc^ , regardless of its proximity to a rich cluster (see also 
Postman fc Geller 198^. 



5<R/R <20 




S(Galaxies Mpc"^) 

Figure 7. The local, projected density distribution of galaxies in 
different radial bins as labelled. Clusters with cr > 800 km s~^ 
are shown as the dotted line, and the mean value is shown in the 
top left corner as Sjjigji. The solid line represents clusters with 
a < 800 km s-^ their mean value is Sio,^. 
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10 1 0.1 



I](Galaxies Mpc ^) 

Figure 8. As Figure ^ but as a function of local projected density. The vertical, dash ed line represents the mean projected density of 
galaxies within the virial radius of the cluster. The solid curves (discussed in Section |4.l[ ) are the expected trends due to the m orphologv- 
density relation of Dressier (1980), assuming the field population is composed of 18% E, 23% SO and 59% spiral galaxies hitmord 
et al. igQ. 




10 1 0.1 



E(Galaxies Mpc ^) 

Figure 9. As Figure 0, but restricted to galaxies beyond 2Rv. 
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4 DISCUSSION 

4.1 Comparison with the morphology-density 
relation 

We have shown that the dependence of star formation rate 
on local galaxy density is independent of cluster ve loci ty dis- 
persion and thus, presumably, mass (see Section 3.3). In a 
recent photometric stu dy based on Hubble Space Telescope 
imaging of 17 clusters, Balogh et al. (2002) found some evi- 
dence that the morphology-density relation does depend on 
cluster X-ray luminosity, which is likely to be a better tracer 
of mass than velocity dispersion; at a given local density, 
low-mass clusters have more disk-dominated galaxies than 
high-mass clusters. Furthermore, they showed that this is 
most likely due to a difference in the population of galaxy 
bulges; the disk luminosity function at a fixed local den- 
sity does not depend on cluster mass. Since star formation 
is generally limited to the galaxy disk, our results are con- 
sistent with this picture. The luminosity of a disk, and its 
star-forming activity, depend only on galaxy density, while 
the luminosity of the bulge component has an additional, 
small dependence on the mass of the embedding structure. 

It would be of great interest to compare the depen- 
dence of /i* on density with the similar density- dependence 
of morphology, to determine the degree to which the two 
correlations are independent. In particular, any difference 
between the two shows that cluster galaxies differ from their 
morphological counterparts in the field, which supports the 
hypo thesis that they hav e undergone a physical transforma- 
tion (Balogh et al. 1998). However, we note that this test is 
not conclusive; if the star formation rate of a spiral galaxy is 
reduced gradually, on timescales similar to that for morpho- 
logical change, the correlation between morphology and star 
formation rate may be retained, despite the transformation. 

Unfortunately, morphological classifications are not yet 
available for our sample. However, we can use the local 



the 41% used here and elsewhere (Whitmore et al. 199J 



Dressier et al. 1997). This high value for the early- type frac- 



morphology-density relation computed by Dressier (1980 ), 
assuming that it is universal. The luminosity limit of our 
sample {Mi, = —19) is similar to that of Dressier, Mi, ~ 
— 19.2, after accounting for the difference in cosmology and 
making the transformation Mh = My + Q.72( B — V) , assum- 
ing a n average galaxy colour B — V = 0.8 (Fukugita et al 



1995). Thus, our density measurements should be compara- 



ble. We will assume that the field galaxy sample is composed 
of 18% E, 23% SO an d 59% spiral and irr egular galaxies 



( [Whitmore et al. 1993| ; [Dressier et al. 1997| ). We therefore 
divide the field galaxy n* distribution (Fig into three 
populations, identifying the lowest 18% of /i* values with 
the E population, the next 23% with the SOs, and the re- 
mainder with spirals. It is then straightforward to recom- 
pute the statistics shown in Fig. ^ for any morphological 
mix. We show the expected /^'-density relation computed in 
this way, assuming Dressler's morphology-density relation, 
as the solid curves in Fig. ^ Two things are immediately 
clear. First, at S = 1 Mpc~^, the lowest density point in 
Dressler's study, the cluster morphological mix is close to 
that adopted for the field, so the predicted curve is in good 
agreement with our measurements. Note that this is depen- 
dent on an accurate determination of the early-type fraction 
in the field, estimates of wh ich have i ncreased from the 20% 
adop ted by Dressie r (1980|), to 30% (Sandage & Tammann 



tion is a consequence of the bright luminosity limit, and is 
consistent with that derived from type-dependent luminos- 
ity functions of [Marzke et al. (1994 40% at M* ). The second 
point is that the predicted ^'-density correlation appears to 
be shallower than the observed relation. This suggests that 
the morphology-density relation may be distinct from the 
star formation-density relation. In making this comparison 
we have made the extreme assumption that the lowest values 
of /i* are associated with elliptical galaxies, and the highest 
values with spiral galaxies. Any dispersion in the natural 
morphology-^* relation will serve to further fiatten the pre- 
dicted ^* -density relation and increase the discrepancy with 
th e data. On the other hand, there is an important caveat, 
as Dressier (198C) did not subdivide the late- type morphol- 
ogy class, and Sa galaxies are known to ha ve much less cur- 
rent star formation than irregular galaxies (Kennicutt 1992; 
Jansen et al. 2000). If the fraction of Sa galaxies relative 



to later types increases with density, this will steepen the 
curves in Fig. H. 



4.2 Possible mechanisms: comparison with 
theoretical models 

These results show conclusively that suppressed star forma- 
tion is not limited to the cores of rich clusters, but is found 
in any environment in which the local projected galaxy den- 
sity exceeds one galaxy brighter than Mb = — 19 per Mpc~'^ 
This is in approximate agreement with the results of 



Ko- 



dama et al. (2001), though a direct comparison is not possi- 



ble because that survey probes much deeper down the lumi- 
nosity function, so the local projected galaxy densities are 
higher in the same environments. Whatever mechanism is 
responsible for terminating star formation in galaxies, then, 
is not particular to the cores of rich clusters, but is associated 
with dense groups in the cluster infall regions as well. This 
means that ram pressure stripping of galaxy disks cannot be 
completely responsible for the correlation of star formation 
with local density, since thi s is only expected to take place 
in the cores of rich clusters ( |Gunn fc Gott 1972 Fujita 2001; 
Quilis et al. 2000[ ). 

Most hierarchical models of galaxy formation do not in- 
clude a calculation of ram-pressure stripping of the cold, disk 
gas, nor of other p hysical processes like galaxy harassment 
(Moore et al. 1999) which might play a role in dense environ- 
ments. The only environmental effect on star formation rate 
in these models - apart from a possible difference in merg- 
ing history - is related to the hot, halo gas hypothesised to 
surround every isolated galaxy. It is assumed that galaxies 
maintain the supply of cold gas - fuel for star formation - 
via continuous cooling from a hot, d iffuse gas halo associ 



ated with the dark matter po t ential ( Bomerville fc Primack 
199S|; [Kauffmann et al. 1999| ; pole et al. 2000(). In haloes 



198l[ ), adopted by [Postman fc Geller (1984 ), and finally to 



with more than one galaxy, this hot gas is only associated 
with the central galaxy; satellite galaxies are assumed to 
lose their supply of fresh fuel through ram pressure stripping 
and tidal effects (though these are not directly modelled). 
In these models, therefore, star formation rates begin to de- 
cline for any satellite galaxy, whether in a poor group or a 
rich cluster. 

These models are able to reproduce radial gradients 
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in star formation within the virial ra dius of clusters to a 
remarkably high degree of a ccuracy ( Diaferi oet_aL_200l| 
Okamoto & Nagashima 2001). In particular, piaferio'et'^r 



(SoOl^JredicF that the mean star formation rate should be 
equivalent to the field value beyond 2Rv, in physical ( i.e. 



not projected) space. The model of Balogh et al. (2000) is 
a greatly simplified version of this more complete model, as 
the properties of the field galaxy population are not mod- 
elled directly, but are taken empirically from observations 
of the z ~ 0.3 field. The advantage is that the effects of the 
halo-stripping can be seen directly, since that is the only 
physical process (apart from gravity) which is accounted 
for. In Fig. we show the predictions of this model, for the 
mean star formation rate relative to the field, as a function of 
local projected density. The simulations on which the model 
is based were kindly provided by Julio Navarro. Here, local 
density is defined as the projected surface mass density, com- 
puted by finding the radius encompassing the ten nearest (in 
projection) particles in th e simulations. The model is the 
"group" model in Fig. 1 of Balogh et al. (2000); galaxies are 
assumed to lose their reservoir of hot gas when they are asso- 
ciated with a group with circular velocity Vc > 600 km s^^. 
While a direct comparison with the data is not possible, 
since these simulations only provide the dark matter density, 
a comparison relative to the mean surface density within 
R-u should be fair if mass traces light. First we note that 
the approximately power-law dependence on local projected 
density has a similar slope in the data and the model; the 
mean star formation rate decreases by a factor of ~ 3 for 
every factor 10 increase in surface density. Secondly, in the 
model the correlation flattens out at surface densities ~ 1/7 
that of the mean projected density within _R„. Although this 
threshold is a factor ~ 2 lower than seen in the data, given 
the crudity of the model, we consider the agreement reassur- 
ing. Unfortunately, the simulations used in this model did 
not include a large enough volume to probe beyond a few 
Ru- Thjis, the low density regions in the simulations are not 
drawn from the same regions m space as the low density 
regions in the observations, most of which are found well 
beyond 2_R„. 

Thus, models in which halo-stripping is the only direct 
environmental-influence on the galaxy star formation rate 
provide a reasonably good match to the data. This is espe- 
cially remarkable given that the stripping is not even directly 
modelled; it is simply assumed that every satellite galaxy has 
no reservoir of hot gas, immediately after it merges with a 
larger halo. Improvement in this respect alone may well im- 
prove the models' success in the lower density regions, far 
from the cluster core. 



4.3 Consequences on the evolution of the global 
star formation rate 

What mechanisms are responsible for driving the stron g ob- 
served evolution of the global star formation rate (e.g. Lilly 



et al. 1996)? One possibility is that the decline in star for- 



mation activity is related to physics internal to individual 
galaxies — for example, consumption of a limited gas sup- 
ply, or a time-dependent cooling rate — regardless of their 
environment. On the other hand, some of the decline is likely 
to be tied to the hierarchical growth of structure; as time 




E (1012 Mpc-2) 



Figure 10. The "group" model of Balogh et al. (2000), in which 
galaxies in groups with circular velocity Vc > 600 km s~l^ have 
their hot gas haloes stripped, so no further cooling is permitted. 
The density is the local projected mass density, computed from 
the area enclosing the nearest ten particles in projection. Small, 
open points are results from a single projection of each of six 
model clusters; the large solid point is the mean, and the error 
bar is the standard deviation of the 18 realizations of the model. 
The dashed line shows the mean density of particles within the 
virial radius. 



progresses, more and more galaxies are locked up in clusters 
where, perhaps, star formation is directly inhibited 



According to exte nded Press-Schechter theory (Bower 



1991 



than 10^** Mq, approximately the limit of our cluster sample, 
is only 11% at the present day, and negligible by 2 = 1. Thus 
it is not immediately obvious that the lower star formation 
rates in these systems can have any effect on the global aver- 
age. However, we have shown that lower star formation rates 
are seen in environments with densities --^ 0.3 times lower 
than the mean cluster density, regardless of their proximity 
to the cluster. This density corresponds approximately to 
the density at the virial radius; by definition, if mass traces 
light then any virialised structure will have a mean density 
which exceeds this threshold. Since our density estimate is 
based on the tenth nearest galaxy brighter than Mt = —19, 
we cannot be sure how our results apply to systems with 
fewer than ten such galaxies. A virialised system with more 
than ten galaxies brighter than this limit is expected to have 
a total gravitational mass M ^ 



Bond et al. 1991 ) the fraction of mass in haloes greater 



mass-to- light ratio of 100 (e.g. Girardi et al. 2002). In con- 



10^^ Mq, assumin g a total 



trast with the more massive clusters, these haloes account 
for ~ 35% of the mass in the present day Universe, and 
contribute significantly to the global average star formation 
rate. Furthermore, at z = 1 only about 10% of the mass was 
in such environments; the rapid growth to z = on these 
mass scales may well be able to explain the rapid evolution 
in the global star formation rate. The hypothesis that the 
growth of structure is largely responsible for the observed 
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decl ine in star formation with cosmic time (e.g. Lilly et al 
1996()lJierefore becomes much more attractive. 



5 CONCLUSIONS 

We have presented a study of seventeen known galaxy clus- 
ters, using redshifts and Ha equivalent widths measured 
from 2dFGRS spectra. We have used this to trace the de- 
pendence of relative star formation rates as a function of 
radius and local density. We conclude the following: 

1. The distribution of star formation rates is correlated 
with both distance from the cluster centre, and with local 
projected density. The distribution becomes equivalent to 
that of the global average for radii ^ 3i?„, and local pro- 
jected densities ^ 0.3 times that of the mean cluster virial- 
ized region. These results are in good agreement wi th pre- 
liminary re sults from the Sloan Digital Sky Survey (Somez 



et al. 2002) 



2. The correlation between star formation rate and local 
projected density holds for galaxies more than two virial 
radii from the cluster centre. Thus, star formation rates de- 
pend primarily on the local density, regardless of their prox- 
imity to a rich cluster. 

3. This means that galaxy transformation is not primarily 
driven by processes like ram pressure stripping, which only 
operate in the most extreme environments, but by processes 
which are effective in lower density, group environments. 

4. The dependence of star formation rate on density is the 
same for clusters with a > 800 km s~^ and for clusters with 
(J < 800 km s~^, which implies that the star formation rate 
is insensitive to the global, large-scale structure in which the 
galaxy is embedded. 

5. The correlation between star formation and density 
that is predict ed from the morphology-density relation of 
Dressier (198C) is less steep than observed. This provides 
conditional support for the view that the correlations with 
density are due to physical transformation of galaxies in 
dense regions, and that morphological change occurs on a 
different timescale from changes to the current star forma- 
tion rate. However, it may also be explained by a lower frac- 
tion of early type spiral galaxies, relative to late types. 
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